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Abstract 

We revisit Heisenberg’s model for nucleon-nucleon scattering which admits 
a saturation of the Froissart bound. We examine its uniqueness, and find 
that up to certain natural generalizations, it is the only action that satu¬ 
rates the bound. We find that we can extract also sub-leading behaviour 
for (Jtot(-s) from it, though that requires a knowledge of the wavefunction so¬ 
lution that is hard to obtain, and a black-disk model allows the calculation 

of (Telasticis) aS Well. 

The wavefunction solution is analyzed perturbatively, and its source is in¬ 
terpreted. Generalizations to several mesons, addition of vector mesons, 
and curved space regimes are also found. We discuss the relations between 
Heisenberg’s model and holographic models that are dual to QCD-like the¬ 
ories. 
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1 Introduction 

In quantum field theory unitarity constraints the asymptotic dependence of the total 
cross section of any scattering process to be bounded by the well known Froissart 
bound [1, 2], 

atot(s) <C'ln2-; C<^, (1.1) 

So 

where s is Mandelstam’s dynamical variable and m is the mass of the lightest particle 
that can be exchanged by the scattering projectiles. In the case of QCD m = 
is the pion mass, and the bound is supposed to be saturated in the s —)■ cx) limit. 
However, being that the saturation of the bound is governed by nonperturbative, 
IR, physics, attempts to describe the saturation of the bound in QCD have not 
been successful. Strikingly, nine years before the discovery of the bound, and in fact 
even before the birth of QCD, Heisenberg proposed a simple effective model for the 
maximal behaviour of atot(>s), in terms of a DBI action for the pion field, that gives 
an almost saturation of the Froissart bound in the case of QCD [3]. 

The paper of Heisenberg [3] includes two revolutionary ideas: (i) Extraction of 
the dependence of the total cross section on Mandelstam s variable from the average 
energy per pion determined from the classical energy density of the scalar field, (ii) 
Describing the dynamics of the scalar field using the DBI action. The first idea is 
obviously very different from the way one usually determines cross section in pertur¬ 
bation theory. Instead of computing Feynman diagrams of scattering amplitudes and 
then from the amplitudes determining the cross section, Heisenberg’s proposition is 
to derive the cross section in a very simple manner from the following relation 

iko) = ^ atot = TTbl,. = ^ log^ , (1.2) 

mi (fco)^ 

where {ko) is the energy per pion, is its mass, hmax is the maximal impact param¬ 
eter for which there is still an interaction between the two nucleon projectiles. The 
assumption of the model is that there in an “effective action” for the scalar field that 
mediates the interaction from which one can compute (/cq) = f where £ is the total 
energy and n is the number of the pions. Thus, the dependence of atot on s follows 
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from the dependence of (/cq) on s. A physical system for which (fco) does not depend 
on s saturates Froissart’s bound. The second original idea is to use a non-standard 
action to describe the dynamics of the pion held. In his paper Heisenberg found that 
using an action of the scalar that is based on ordinary kinetic term and regardless 
of what is its potential cannot saturate the Froissart bound. In fact, it will yield a 
constant cross section. However, using a DBI action yields (fco) ~ log This mild 
dependence of (ko) on s means that the total cross section of the model is close to 
that of the bound. 

Experimental data of the total cross section of proton-proton (and proton- an¬ 
tiproton) collisions is well established on a very wide range of energies starting from 
sub GEV energies and all the way to ^/s = 7 in the TOTEM experiment in the LHC 
and in fact even higher up to ^/s = 57 TeV from cosmic rays observation. Figure 
(1) shows the data points together with a ht based on the (1.2) but with a mass 
m ~ IGev and not the pion mass [4]. Thus, regardless of the Froissart bound, one 
would like to have a theoretical model that resembles the behavior of (1.2) since it 
seems to ht the experimental data quite well. Needless to say that there is no direct 
derivation from QCD that can reproduce such a ht. 



Figure 1: The total cross section as a function of for pp and pp scattering 

The goals of this paper are fourfold: (i) Since part of the paper of [3] is written in a 
concise form, we elaborate the discussion, performing several additional calculations 
and provide some further evidence for the claims of the paper. In particular we 
analyze the pion held including the passage from the 1-1-1 dimensional solution to 
a full four dimensional one. The ratio of elastic to total cross section is derived 
using a black disk model, (ii) We examine the uniqueness of the DBI action as 
an action that can (almost) saturate the bound. We prove that only an action 
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with an infinite tower of higher powers of the derivative term, as the DBI action 
admits, can do the job. (iii) We propose and analyze several generalizations of 
Heisenberg’s model. We add a general potential instead of only a mass term, and 
we analyze a sigma model with several scalars. We examine the “highly effective 
action” of [5]for the case of single scalar in AdS^. For that case the ordinary kinetic 
term in the square root Lagrangian density undergoes the following transformation 
df^(j)d^(j) —)■ Upon considering an power rather than 0^, we show 

that only for the range n G (0, 2) we obtain saturation of the Froissart bound, 
(iv) The last goal has been to relate Heisenberg’s model to the DBI action used in 
gauge/gravity duality and furthermore to two different holographic approaches to 
the nucleon-nucleon scattering. 

One approach to the latter is based on a simple effective model for QCD scattering 
at high energies has been developed, the Polchinski-Strassler model [6] in terms of 
a metric in a cut-off HdS's background, dual to glueball fields, and the fluctuation 
of an IR brane (IR cut-off), dual to a pion field (the model was extended to the 
phenomenological ’’hard-wall” model with the addition -by hand- of extra fields in 
the bulk in [7]). Based on earlier work in [8], in [9] it was shown that the saturation of 
the Froissart bound can arise through scattering of gravitational shockwaves located 
on (or close to) the IR brane, with formation of black hole on the IR brane. Then 
in [10] it was shown that one can map exactly the description of the saturation of 
the bound in the dual, through gravitational shockwave collision with black hole 
formation, with the saturation of the bound in the Heisenberg model, through pion 
field shockwave collisions. This picture was used further in [11, 12, 13] (see [14] for a 
review) to describe the sQGP fireball obtained in heavy ion collisions as the object 
dual to the black hole formed on the IR brane. 

In the second approach we relate Heisenberg’s DBI action to the DBI that de¬ 
scribes the fluctuation of the flavor branes in confining gravitational background 
and in particular in the generalized Sakai-Sugimoto model [15] [16]. We show that 
what sources the scalar field in that model is the flavor instanton density [17], that 
corresponds to the proton density. 

The paper is organized as follows. In the next section we review the model of 
[3] and elaborate on certain issues by performing additional calculations. In section 
3 we discuss the uniqueness of the DBI action used in [3] for the saturation of the 
Froissart bound. We show that to have a solution of the form ^(s) = A^/s, which 
is what is needed to saturate the Froissart bound, one cannot use an action with a 
finite series of higher derivatives and only the infinite series that follows from the DBI 
action does the job (though we were not able to show that some other action with an 
infinite number of higher derivative terms couldn’t do the job also). In section 4 we 


4 



consider various generalization of the action used in [3]. First we add a potential term 
in addition to the mass term. We then analyze a DBI sigma model for several scalar 
helds. For the simple case of replacing one kinetic term in the square root with a 
sum of kinetic terms for several scalars the behaviour is similar to the original model. 
Next we discuss the case of a DBI action associated with the AdS^ x case. We 
show that for the case of a single scalar held performing the determinant in the DBI 
action yields a close cousin of action used in [3]. Another generalization discussed is 
the DBI for vector mesons. Assuming here again dependence only on the coordinate 
s dehned in (2.3) we show that the behaviour of the vector mesons is similar to that of 
the scalar meson. In section 5 we discuss the “wavefunction” of the pion. Firstly we 
elevate the solution ^(s) to 0(s,r). We argue that for Heisenberg’s solution we have 
a shock-like behaviour, where T++ blows up at = 0, even though we don’t have 
a behaviour. Next we consider a perturbative expansion around r = 0 and 

then an pertubative solution around the asymptotics r —)■ cx). Section 6 is devoted 
to analyzing the sources of the pion held. We hrst consider sources for the nonlinear 
Born Infeld theory of electrodynamics. We then in subsection 6.2 discuss in a similar 
manner the source of a scalar DBI theory. Section 7 deals with the original question 
behind this paper, namely, the cross section of the nucleon-nucleon scattering process. 
We discuss corrections away from the Froissart bound. We then describe the model 
of the black disk and the corresponding ratio between the elastic and total cross 
sections. Section 8 is devoted to an examination of the relation between Heisenberg’ 
model and the holographic description of nucleon-nucleon scattering. We show that 
the nucleon-nucleon scattering process takes the form of the scattering instantons of 
the flavored gauge helds that reside on the havor branes. We end this paper in a 
section of summary and open questions. 


2 Heisenberg’s model- A review and elaborations 


In this section, we hrst review the work of Heisenberg in modern language, and we 
perform some additional computations that clarify some aspects of the model. 

With a remarkable insight, Heisenberg considered a nonlinear higher derivative 
action for the pion, the DBI action with a mass term inside the square root. 


C 


7-4 


1 - y7+77777777Fj 


( 2 . 1 ) 


The reasoning is that in the high energy limit, many pions (lowest mass particles) will 
be created, so we need to consider a pion held as the ehective one, but this process 
is both nonperturbative and high energy, hence one needs a nonlinear action. As we 
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will soon see, a polynomial interaction does not have the required properties, so the 
DBI action is the natural one to consider. 

In the high energy limit, colliding hadrons will look like pancakes due to Lorentz 
contraction, but moreover we need to consider them as just sources for the pion held 
surrounding them, that will also get Lorentz contracted and look like a shockwave. 
Therefore the process considered in the asymptotic regime is a collision of pion held 
shockwaves with the action (2.1). We look for (classical) shockwave solutions to the 
action (2.1). The equations of motion are 


-□0 + + 1 ‘ 


[{d^d^(p){df,(j))d^4> + (( 9 ^ 0 ) 

1 + /"‘[(( 9 ^ 0)2 + 771 ^ 02 ] 


= 0 . 


( 2 . 2 ) 


The crucial simplihcation that allowed Heisenberg to do exact calculations is to 
consider that for a shockwave solution, only the physics near the shock is relevant, and 
by focusing near that, we can ignore the dependence on the 2 transverse dimensions 
y,z (with r = \/y‘^ + z'^), and consider the 1+1 dimensional problem for time t and 
longitudinal direction x (along the direction of propagation). 

Then from Lorentz invariance, he considers only solutions that depend on 

s = t^-x‘^. ( 2 . 3 ) 

(Note that from now on, we will use s to denote this variable only, and not the 
Mandelstam invariant, which will be called s) This requires some explanation. The 
hrst point is that 0 = 0(s) is boost invariant for boosts in x: Under a boost, we have 
x~^ — )■ e^x~^,x~ — )■ e~^x~. But why do we need a boost invariant solution? The fact 
that 0 is a scalar means that (l)'{x'~^,x'~) = (l){x~^,x~), where = t ±x. We could 
say that we hnd the solution 4>{x^, x~) in a reference system and then dehne the one 
in another reference system by (j)'{x'^, x'~) = 0(a;+, x~), so any solution would work. 

However, the essential point is that we use the ultra-relativistic approximation, 
in which even though the pion is massive, we consider that the source moves on 
a lightcone, = 0 or = 0. As a result, we impose that (j){x'^ = 0) = 0 or 
0(a;“ = 0) = 0. This in turn implies a power-law behaviour near the lightcone, i.e. 
(for x~ = 0), 0 oc {x~y, g > 0 for ~ 0. But if we have an arbitrary dependence 
on then in the boosted system 0' would have a power of in front, unless we 
have the same power law for x'^ , i.e. unless 0(a;+,a;“) = (l){x^x~) = 0(s). 

For 0 = 0(s), we have 
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the DBI action becomes 


-4 



1 + /M -4s ( 

L \ 

v ^ 


C = I 

and its equation of motion becomes 


dcj) 

ds 


i2/A2 


(2.5) 


4if-(s^)+„2A + 


ds V ds 


l + /4[-4s(f) +m202] L 


d(j) d'^cj) 

ds 2 ds'^ 


= 0 . ( 2 . 6 ) 


However, by multiplying with the denominator (assuming that it does not vanish), 
canceling and rewriting the terms we are led to the form 


d f d(j) 


ds V ds 


4 ^ = 8sl^ ^ 




ds J 1 + l^rn^cj)^ 


(2.7) 


When m = 0, one can hnd an exact solution depending on an arbitrary parameter 


a. 


and 0 = 0 for s < 0. 

When m 7^ 0, one can find a perturbative solution at small s, 

fs 

0 = ^(1 + a s rri^ + ...) , 0 < s < 1/m^ , 

r 

and 0 = 0 for s < 0, as well as a solution at large s, 

(j) ~ cos{m\/s + 5) , s 3> 1/m^. 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 


At this point, Heisenberg notes that for the model to be reasonable, we need that 
{d^4>y to be a hnite constant at the position of the shock, s = 0, since we need the 
nonlinearities to play a role there. But for the free KG equation, the result is inhnite, 
which is also unphysical.^ Thus the only possibility to correctly describe the shock 
at s = 0 is to have (i9^0)^ a hnite constant, which leads to 0 ~ A^/s for s —)■ 0. This, 
as we will show below, is incompatible with an action with a canonical kinetic term 
and a polynomial potential. 

^The KG equation for (j) = 4’{s) is just d/ds{s d<j)/ds) = 0, with the solution (j) = Alog(s/so), 
which means (d^b)^ = —4,s{d4>/ds)'^ = —AA^js —)■ oo. 
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2.1 From the pion field to the nucleon-nucleon cross section 


The energy (Hamiltonian) density of the pion held is 


H = 7r0 — £ 


a/1 + H [((9^0)2 + 


and similarly the momentum density is 


( 2 . 11 ) 


V = 7rV0 


_0V0_ 

a/I + /"^[(i9^0)^ + 


( 2 . 12 ) 


Both densities have a denominator which is the square root term of the Lagrangian 
density. In the massless case, substituting the solution (2.8) into (2.11) we hud 
that the energy density diverges at s = 0 due to the denominator going to zero 
as a^/s/(2l‘^). Similarly, in the massive case, substituting the solution (2.9) into 
(2.11) we hnd the same divergence due to the denominator going to zero at s = 0 as 
m^/s{l — 6a). 

Now following [3] we assume that one can introduce a small perturbation so that 
the denominator can be taken as a non-vanishing constant. In this case we can use 
the standard method of Fourier transforming (2.9) over x to k, as 


= I ^ — x‘^{l -I- am^it^ — x^) d-...) 


(2.13) 


which for a = 0 (only the leading term) gives 


(j){k, t) ^ I 


2\k\ 


(J.(|*:||(|)+*Hi(|k||t|)), 


(2.14) 


where Ji is a Bessel function and Hi is a Struve function. When expanded at large 
k, we obtain 


4, - (i + . (2.15) 

Note that the non-oscillatory part of (j) is not a radiative piece, hence is dropped. 

However, as discussed above, in reality the shockwave should have a hnite thick¬ 
ness in yds of the order of the Lorentz contracted 1/m, i.e. v^om = V^min ~ 
v‘^/m, which means that at sufficiently large t, (f){k^t) should be cut off at 
korn = the relativistic mass of the pion. With the assumption of a 

constant denominator we get 


dk 


(X k‘^4>{kY 


const. 


(2.16) 


rsj 


k 



where in the last equality we have substituted (2.15). But this is valid only for 
k ^ ^Om- 

Finally, the momentum k is identihed with the momentum of a pion ko, and 
moreover the classical held close to the shock is identihed with the classical limit 
of the held of radiated pions in a hadron collision. Thus the radiated energy £ 
(identihed through canonical quantization with the pion held energy E) per unit 
frequency of radiated pions is given by (denoting the constant by B) 


d£ _ B 

dko ko ’ 


m < k < kom- 


(2.17) 


This integrates to 

£ = B\n— = B\n-f , (2.18) 

m 

and leads to a relation for the number of pions emitted for a given energy, since 
dE = kodn, giving 

(ItIj B 

— = 72 1 m<ko<kom , (2.19) 

a/Co /Cq 


which integrates to 



Then the average emitted energy per pion is 


( 2 . 20 ) 


{ko) 


£ \n{kom/m) In 7 

— = m - - — = m - Y — ■'^Iny , 

n 1 - m kom 1-7 


( 2 . 21 ) 


which is approximately constant (only logarithmic dependence on the energy). 

The last step in the Heisenberg model is to assume that the emitted energy is 
proportional to the total energy of the system, (here s is the Mandelstam vari¬ 
able), with the constant of proportionality (ratio of emitted energy) being approx¬ 
imately given by the pion wavefunction overlap. Since at large transverse distance 
r (= + z"^), the wavefunction is small 0/ 1, thus it satishes the free massive 

KG equation, with solution 0(r) ~ the wavefunction overlap is ~ 6“™^, where 

b is the impact parameter, i.e. transverse separation between the colliding hadrons 
at the impact point a; = 0. Then we have approximately 




( 2 . 22 ) 


The maximum impact parameter for which we have interaction, femax, arises when the 
emitted energy equals the average emitted energy per pion (ko), so that it corresponds 
to emitting just one pion. Then we have 

= {ko)^br.,^ = -\n^^ 

m {ko) 
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TT 9 \ S , > 

atot = 2.23 

{ko} 

We see then that the saturation of the Froissart bound arises only if {ko) is approxi¬ 
mately constant as a function of energy. 

Next we would like to compare this result with what one gets for an “ordinary 
held theory” with a canonical kinetic term and a polynomial potential of the form 

C =- X(j)^ , (2.24) 

The corresponding equation of motion resulting from it for the 0 = (j){s) ansatz, 

^ = 0 , (2.25) 


do not have 0 ~ A^/s as a solution, since the hrst bracket is divergent, as it equals 
—A/2y/s, and the other terms give zero, as they are positive powers of s. In fact, 
we can see that the only way to satisfy the equation of motion at leading order in 
s with a canonical kinetic term plus a potential is for a potential that includes the 
logarithmic term Aln(0/(;/)o), since then in the equation of motion we have A/^, and 
we can solve the equation with A^ = 2A. But it is unclear how such a term could 
arise in the potential (especially since it is unbounded from below at 0 = 0). 

Now let’s study for this class of theories the energy per emitted pion. In a way 
similar to the one described above we can prove that 


dS 

dko 

(ko) 


S 

_ r 

n 


m<ko< kom] ^ 


dn 

dko 


B 


kom 1 V ^ 

—r— = my--oc- 

In^ Iny InVs 


m <ko<kom^ 


(2.26) 


That means that we don’t get the saturation of the Froissart bound, but rather we 
get a constant o'tot(V^)- 

In fact, we can check that the saturation of the bound is obtained only for 
dS/dko oc l/Zcg, with n > 1, whereas for actions with polynomial potentials this 
is not satisfied. In fact, as we saw, dS/dko oc l/Zco was obtained from the behaviour 
0 oc ^/s of the held near s = 0, which was due to the DBI form of the action. 

In conclusion, we have two physical ways to restrict the form of the action. As 
Heisenberg argued, we need (0^0)^ to be a hnite constant in order to describe the 
correct physics, which restricts to 0(s) oc ^/s, arising only in DBI. On the other 
hand, if we are to be able to saturate the Froissart bound (which should happen, 
as Froissart argued), we again need dS/dko oc l/ko, which again requires the DBI 
action. 
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3 Uniqueness of the Heisenberg model action. 


In his paper [3], Heisenberg shows that an action based on ordinary kinetic term with 
any kind of a potential term does not saturate the Froissart bound. We now seek to 
check how unique is the choice of Heisenberg of having the DBI action as the action 
of the pion held. We have seen that we need an action with higher derivatives. The 
question is then can we have other higher derivative actions? In particular we want 
to examine whether one needs an inhnite series of any power of the derivative term 
or it is enough to have certain hnite series. And furthermore if one needs an inhnite 
series is the DBI action used by Heisenberg unique? 

We now examine this question, by considering Lagrangeans of the type £(0, X), 
where X = {8^4)^. 

DBI truncated to first term. 

We will start by truncating the DBI Lagrangean from the previous section (with 
promoted to 2V for more generality) to the hrst interaction term, i.e. 

= + + ( 3 . 1 ) 

where V{(j)) = ^(0) ~ but for generality we will consider arbitrary V. 

Then the equation of motion is 

( 3 . 2 ) 

and for 0 = 0(s) we get the equation of motion 




■/V(0)) + W(0) + 8s/^ ( ^ 


cPcj) 




= 0. (3.3) 


We want to see whether a solution of the type 0 = A^/s near s = 0 is possible. 
First we note that in this case, the terms with V and V are irrelevant (they are 
subleading), so we will drop them for simplicity (we can add them for free at the 
end). Then, substituting, we get the equation of motion for the leading term 


A 





0 , 


(3.4) 


so we see that for a real scalar held (as we want), when A^ > 0, there is no solution. 

Generalization with first derivative interaction 
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Next, we drop the irrelevant V terms and generalize by writing an arbitrary 
coefficient for the interaction term, 

£ = + C[(a,0)T . (3.5) 

giving the equation of motion 

-□0 + 4C(d^(/>fn(/> + 8C(d^(/>)(d,(/>)(d^d,(/>) = 0 , (3.6) 


and on 0 = 0(s), we get 

ds ^ds^ 


d ( d(j) 


ds V ds 




d(j)\ d f d(j) 


ds J ds \ ds 


+ 64Cs ( ^ 
ds 


0. (3.7) 


Substituting the ansatz 0 ~ A^/s, we get 


+ (3.8) 

Note that we have kept the last term, which is zero, for reasons to be explained later. 

Generalization to arbitrary powers 

Next we consider on top of the previous, an arbitrary n-th order interaction, 

£ = + C2l(9„0)y + C„|(a,..^)y , (3.9) 

with equation of motion 


-□0 + 4C'2(9^0)'n0 


+8G2 {du(i>) {d^dy(t)) 


1 + 

n{n — 1)G„ 


1 + 


2G2 


[(•9^0) 


2in —2 


= 0 , 


(3.10) 


and on the ansatz 0 = ^(s) = A^/s, we have 


/A A 


+16^271" 


\2y^ 2^/s 


1 + 


n{n l)Cn / _^2^r^-2 


2Go 


= 0 . 


(3.11) 


We can hnally generalize to a sum of arbitrary powers. 




(3.12) 


n>2 
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with equation of motion 

-□0 + + ((9^0)(<9^0)((9^(9^0) ^4n(n - 1)[((9^0)^]”“^ = 0. 

n>2 n>2 

(3.13) 

On the solution 0 = 0(s) = Ay/s, we get 


A 

7^ 


n + '^2nCn(-irA^^"-^n 

V n>2 / 


f A A 


^8n(n-l)C'n(-l)M2(”-^) = 0. 

n>2 

(3.14) 

For the DBI action, the coefficients, coming from the expansion of 


.(1 + 5 _ l/2(l/2-l) ..(l/2 - n + l) ^„ 

2 1 * 2 * ... * 


(3.15) 


give therefore sgn{Cn) = (—1)”, meaning that the coefficients inside the two brackets 
in (3.14) are all positive. Moreover, from the arguments in [18], the signs coming from 
the DBI action are the ones needed for causality and locality of an action (note that 
the metric convention in there is mostly minus, so all coefficients there are positive). 

That means that for a general action, at any hnite order in the terms, 0 = Ay/s 
is not a solution. 

But then the question is, how is it possible that the DBI action has this as a 
solution? To answer that, we look at equation (2.6), which is the equivalent of what 
we have here. On the ansatz 0 = Ay/s, we get from it 


A 2/^2 / A ^ ^ _ 

7s 1 - A‘‘p (Vs ~ ~ ° ’ 


(3.16) 


which at hrst seems not to have a solution, just like our hnite order truncations, but 
looking better we see that A^ = is a solution, since then the second term is 0/0, 
and there is a solution, as seen by going to the form (2.7). The essential fact is the 
existence of the factor 

—^— = 1 + X + x'^ + ... + x'^ + ... (3.17) 

1 — X 


for X = 1, multiplying the {A/2y/s — A/2y/s) = 0 term, but not the nonzero term. 
Thus in the case of the hnite truncation, we have the ratio of the zero and nonzero 
terms being 


E„>2M^-l)C'n(-l)M2(-l) 
l + En>2 2na(-l)M2(-i) 


which goes to inhnity for an inhnite number of terms, allowing the solution. 


13 



In conclusion, the DBI action is the unique one satisfying the physical requirement 
0(s) ~ Ay/s near s = 0 (it could be that there are other derivative actions, with 
an inhnite number of terms, and the same signs as DBI for the coefficients, but it is 
unlikely), however we can add a potential inside or outside the square root without 
modifying the result. 


4 Generalizations 


We hrst consider a simple generalization, instead of just a mass term inside the 
square root, a general potential V, with Lagrangean 


C 




1_^1 + / 4 [( 5 ^ 0)2 + 2 D ( 0 )] 


(4.1) 


Its equation of motion is 


-□0 + a^D(0) + / 


4 [(<9^<9^0)(<9^0)<9^0 + {d^(f))^d^V (0)] 
l + /4[(a^0)2 + 2D(0)] 


= 0 


(4.2) 


and for a solution 0 = 0(s), we obtain (after the same manipulations as in the 
Heisenberg case) 


4 


A. 

ds 




(d4>y [f.+v'm 

v*/ 1 + 21‘‘V(4j)' 


(4.3) 


It is easy to check that this equation has again the same small s solution (2.9) for 
a = 0, i.e. the leading term, since the terms with V in the equation of motion 
are actually subleading with respect to the others. This in turn leads to the same 
analysis of Heisenberg, so this generalization is allowed. 

We can also consider adding V outside the square root. 


C 





vw, 


(4.4) 


and we can again check that the same thing happens: the solution (j){s) = l~‘^y/s +... 
is still valid, since again the terms with V in the equation of motion are subleading 
on the solution. 


4.1 Several mesons and sigma model 

We can also consider N scalar helds, corresponding to having several scalar mesons, 
0b i = 1, ...,N and for generality consider it in d + 1 dimensions. A generalized DBI 
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model would be 


£ = + 2V{4>^)] . (4.5) 

Then when /i(0*) = (7(0*) = 1, for small helds we keep only the leading term in the 
expansion of the square root, and obtain the usual sigma model with a potential, 

£ 2 « - VW)- (4.6) 

The equations of motion of the action (4.5) are 

d^ih{4>) - l<^+^[gij{4>^){d^4>%d^4>3) - 21/(00] x 

X v^l -/'^+i[^y(0^)((9^00((9^0J') -21/(00] + 

/(0) ^'9‘^‘£/ifc(0)(0/.0^)0^0^ - 00414(0) - (0)0^0^] - dj{(j))gij{(j))df,(j)^ 

_hid+l) _ /( 0 )^*^'( 0 ) 0 ^ 0 ^' _ 

2 1 - /'^+l [^y(0^)(0^00(0^0t) - 21/(00] 

>< (0M[^n(</>)(0M<^*)(0M</>^) - 214(0*)]) = 0. 

(4.7) 

To analyze this case, hrst note that Heisenberg already considered the case of 
several mesons with DBI action, but that meant that the sum was outside the square 
root, 

C = l-^Y^ 1 - y^l + H[(0^0“)2 + m202] . (4.8) 

a *- -* 

That case worked in the same way as for a single meson. We now consider the 
generalization with the sum inside the square root, and a sigma model metric, 

C = l-^ 1 - . 1 + H ^G'ab(0c)(0;.0“)(0^00 + , (4.9) 

\ L j 

where we can replace everywhere the mass terms with a general potential, since as 
we already saw that doesn’t change anything. 

Now in terms of the asymptotic value of the cross section (the Froissart be¬ 
haviour), nothing changes, since the maximum cross section is governed by the pion 
of smallest mass, that has the largest wavefunction at large distances, according to 
the mechanism reviewed below. What does change is the value of the cross section 
at intermediate energies, where now we have cross sections for emissions of various 
scalar mesons. 
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The equations of motion coming from the action (4.9) are 
d 

- I 

ds 

dcj)^ d<p' 


4^ ( ) Tm" 


2 

a'f'a 


1 + /' 


, ^ dcj)^ dfj)^ 2 , 2 


-fmtG. 


ab 


-As- 


ds ds 


r - ssi^Gab 


d<p° d 
ds ds 


ds ds 
d(j)^ d(l)^ 
ds ds 


= 0 . (4.10) 


The simplihcations that occured when Gab = dab do not occur anymore. However, 
the helds 0“ will have in general also the interpretation of some brane coordinates 
in the gravity dual descriptions of section 8. In the Heisenberg case, we had a single 
held, corresponding to a single coordinate transverse to the brane, but in general 
we can have many. Then the origin of coordinate, corresponding to the position 
of the brane, must be a stable point. Around it, we can expand the metric as 
Gab = (^afe + C>(|0|), and write an ansatz for the helds as 

0“ = , (4.11) 

for s —)■ 0. At s = 0, the helds are at 0, i.e. the stable point (the ”IR brane” or IR 
cut-oh of the gravity dual), and with the metric expanded as above, we can check 
that the ansatz is a solution of the equations of motion if 

^(A“)2 = l/l\ (4.12) 


In order to understand the asymptotic cross sections, we consider the behaviour 
of the wavefunctions for large transverse r. The large r behaviour of the cross 
section is governed by the lightest meson, the pion. Indeed, the wavefunctions go 
like 0“ oc 6“™“” at r —)■ cx). Therefore, if we are in the asymptotic regime for the 
pion, 0^ oc which implies also |0“| \(jf\ then from the equations of 

motion we can check that we also have 0“ oc which by the usual Heisenberg 

argument implies that 


6 ~ 


^-maba ^ 


(kp) 

V~s 

(M 


^ . 2 f 

atot ~ ~ ^ In 

mi V {kp) 

/ I — \ 

^ 1 2 
0-a 2^ ^ In 

mi 


(^o) 


^tot 7 


(4.13) 


where a a is the cross section for production of mesons a. Therefore in the asymp¬ 
totic regime, all the Oa should behave like Froissart saturation, with corresponding 
coefficients 7r/m^. 


4.2 AdS case and curved space generalizations; solutions 

In section 8 we will discuss possible relations between Heisenberg model and a holo¬ 
graphic description of the nucleon-nucleon scattering. We have seen in the previous 
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subsection that a natural generalization of Heisenberg’s model includes several scalars 
corresponding to several mesons, a sigma model for them. Here we discuss a partic¬ 
ular example of such a generalization which is a DBI sigma model in AdS spacetime. 
This arises also naturally in the context of gauge/gravity duality. 

The DBI action on the flat worldvolume in d -|- 1 dimensions, i.e. for a Dd-brane, 
takes the form 

Sdbi = TdJ d"+Ve-y- det[d^X^d,X^gijiX) + 2na’F^,] , (4.14) 


where is the Dd-brane tension, 0 is the dilaton, af^ are the world volume coor¬ 
dinates, X* i = are the target space coordinates and gij{X^) is the metric 

on that target space. The d -|- 1 dimensional DBI action describes in particular the 
physics of Dd-branes. The D-brane action in fact also include a CS term, but for our 
purposes, the effect of that will be just to subtract J from the above action. 

Imposing d -|- 1 dimensional Lorentz invariance, switching off the gauge helds, 
writing using the static gauge for / = 0,1, 2, d, and dehning 

the vector 0 = = in* with i = d -|- 1,..., D, the DBI action reduces to 


Sod = / d"+' 


xe 


- det irii„yg{(j)) + /'^+^<9^0*<9,.0^fi'ij(0)) - 1 


(4.15) 


For the special case of the D3-brane moving in an AdS^ x space (the space 
generated by a large number N of other D3-branes), this action is the ’’highly effective 
action” for the A/” = 4 SYM theory on D3-branes written recently in [5], which takes 
the form 


A 


D3 


d^a;0^ 




-det -f 


< 9^0 ■ < 9^0 


- 1 


(4.16) 


where now the dilaton is a constant, the target space coordinates are 0 = n* with 
I = 4,..., 9, and the metric on the target space was taken to be 


ds^ 




(f^gijdx^ dx"^ + 


—d0 ■ d0 


(4.17) 


For the special case of a single scalar with 0 = 0, ^ = 1, using the identity 

- det + g{(j))df,(j)d^4>] = 1 -h g{(j))df,4>d>^4 >, (4.18) 

the action for a D-brane (4.15) reduces to the Heisenberg type action with a g{4>), 
or one-dimensional sigma model action of the type (4.9). 
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For the case of several scalars, but with the metric trivialized around the position 
of the brane, i.e. gij{(t)) — g{(l))Sij, we have 

- det(r/^^ + g{(j))d^(j)^d^(j)^) = <9^,10*)... 

+...] 

= 1 + g{4>)d,4>^dy^ + 

+g^ (0) (<9^10*' <9^10*' ) (<9^2 0*" <9^2 0*' ) 0*' <91^3 ) <^w)l2'^'3 

+... (4.19) 


Then on the solution 0 = 0(s), we have 


2 (<9/,i 0*) (9^2 0^ (9^2 ) 5); 

= ( 9 ^^ 0 *( 9 ^^ 0 *( 9^2 0 -^( 9^2 0 '^ — ( 9^2 0 * ( 9^2 


^IJ-lU2 


= 16s 


2 d(j)'‘ (i0* d(f)^ d4d ^g2 dcf)^ d^d 


= 0 


(4.20) 

ds ds ds ds ds ds ds ds 

and we can easily see that for the higher terms the same happens. Therefore on the 
solution 0 = 0(s), the presence of higher order terms inside the square root in the 
D-brane DBI action (4.15) is not relevant, and we have still a sigma model action 
like (4.9). That means that the Heisenberg analysis is still valid on the case of the 
general DBI D-brane action. 

Shockwaves for D-brane in cnrved space 

Consider the action of a D3-brane moving in AdS^, i.e. the ’’highly effective 
action” of (4.16) for a single scalar 0 and metric g{(j)) = 0“^, 


C = l 


-4 


1 + 




Its equation of motion on the ansatz 0 = 0(s) is 

.0" + 0'-2^(0')2-2^(0')3 = O, 


(4.21) 


(4.22) 


which is a special case of the more general form with an arbitrary metric g{4>), 


s(j)" + 4>' + 2 ~ 2'Sfl'(0)(0O^ — 0 , 

for g{(j)) = 

It is easy to check that an exact solution for this non-linear equation is 

<l>{s) = ^ = ^ 


(4.23) 


'/s y/t‘^ — x'^ 


(4.24) 
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Substituting the solution into the Lagrangian density (4.21) we find that the 
square root vanishes and C = (The same holds for the solution of the massless 
Heisenberg model where g{(j)) = 1). 

In fact one can use this property to find solntions for other target space metrics. 
For g{(j)) = we get 


C = l-^ 


4s^(0)(0')^ = 1 ^ 


0(s) = / ^(/ ^s)2-» 


2 — n 
4 


1 

2 — n 


(4.25) 


Furthermore, for a general metric g{(j)) we find that the solution for 0(s) is 


J d(t)\fg{^ = ^/s. 


(4.26) 


It is easy to check that this solution solves indeed the equation of motion (4.23). 

We now want to see whether we can satnrate the Froissart bonnd for 5f(0) = 0“^, 
with solntion (4.24). 

Its Fonrier transform is 


0(f,/c) = 


-)■ 


^{M\k\\t\)+Lomm 

/¥ (l + 5^) 


(4.27) 


where on the second line we have written the /c —)■ cx) limit, and Lq is a Struve 
function. 


Then for the energy per pion frequency, we get 


dS 

dko 


k‘^\(l){k, t)p ~ k. 


(4.28) 


As we already argued, in this case we do not get a satnration of the Froissart bound. 

Similarly, for g{(j)) = the Fourier transform of the solntion (4.25) gives 

at K —)■ cx). 


(pit, k) — il 2- 


Xi 


n —3 
X \ n-2 


2 — n \ 


_ 4 —n . . j —n . . 1 _ 1 

—)■ / 2-„(^_^j2-n(^2 — n)2-n2 2-n T 


n — 3 


X 


f 1 + 

\k\ 


n — 3 


2(n — 2yi \k\ \t\J 


n — 2^ 
(4.29) 


This leads to 


^ = k‘^\(l)it,k)\^ ~ A;"-2 , 

dko 


which means we obtain satnration of the bound only for n G (0,2). 


(4.30) 


19 



4.3 Introducing vector mesons 


As we saw in (4.14), the DBI action on the worldvolume of a Dd-brane has vector 
helds. In fact, in AdS/QCD approaches with probe branes, like for instance the 
Sakai-Sugimoto model or the models of section 8, these vectors on the gravitational 
side give rise on the dual held theory side to towers of vector meson states. We are 
considering mainly a hat metric for the scalar helds (trivial sigma model), which 
arises as an approximation in the IR of the gravity dual, as we discussed. Then 
we must consider the ehect of the gravity dual metric (that drives the brane to the 
stable point around which the metric is hat) to be to give masses to the helds. 

Therefore we consider the DBI action with a mass for the vector inside the square 
root, 

c = I 

= 1 


where Fab = ^eabcdF^^ . 

For the vector wavefunctions Aa{r), like for the pion held 0(r), we need to give 
some initial data (boundary condition), and then the wavefunction is determined 
from the equation of motion of the above action. 

Let us consider hrst the case with no pions, just vector mesons, i.e. 0 = 0. At 
sufficiently large r we have again the usual free held decay 

Aa{r) = Aae-^^^ , (4.32) 

and again, with the additional assumption that (jy^ the cross section for emission 
of V vector mesons, is obtained when the emitted vector meson energy (which we 
should calculate) equals the average per vector meson emitted energy, i.e. 

g ^V^max 

we obtain 

= tt&Lx- (4-34) 

Of course, the correct calculation would be the one where we have both the pions 
and the vector meson wavefunctions, and then we can calculate ay as above. 


(4.33) 
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The action for only vector mesons with mass My and no pions is 


C = I 


-4 


1 - 


\ 


1 + 


+ 


(4.35) 


We would like to restrict again the dependence of the gauge helds to a dependence 
on {t, x) and furthermore to only s dependence, but now with all the four vector helds 
Aa, with a = 0,1,2,3, since there is no gauge invariance due to the fact that the 
vector mesons are massive ones. 


Substituting Aa = Aa('S) in Fab-, we hnd 




= -8[s 
1 


dA-^ 

ds 


+ s 
1 


dA. 

ds 


dAi dAn 

+ \t^ + x 


ds 


ds 


ds ds ds ds 
so that the action for the ansatz is 


r rjiaft _ ^ abed TA ta _ ^^ 0123 /qt 7 ta oi?i? 

XabX — -e rab-Tcd—-^^ 10^01-^23 ” 0 ^ 02.^13 + 03-^12J 

Q dA2 dAs dA^ dA2 . 

—ot— — X— -h t— — X— — 1=0, 


dAi dA^ 

+ t— -h X 


(4.36) 


C = I 


-4 


1-4/^ 


dA2 

ds 


+ s 
1/2' 


dA^ 

ds 


ds 


ds 


+FM^{Af, + 


(4.37) 


Note that in this action we can consistently truncate Aq = Ai = 0, and then the 
DBI action for A 2 and A^ are the same as for two DBI pions of Heisenberg, for which 
we already saw that we need the full nonlinear DBI action. 


5 The pion wavefunction 

The pion wavefunction should be a solution of the equations of motion coming from 
the pion action. Following Heisenberg, we have considered only the 1+1 dimensional 
case of (f){s) that describes the physics near the shock, at s ~ 0, and the weak held 
case 4>{r), spherically symmetric in the transverse coordinates, so a function of only 
r = a/i/^ + z'^. 

Note that in general, we do not need to have even an ansatz depending on both r 
and s, i.e. (j){s,r), but rather depending independently on all 4 coordinates, however 
considering (j){s, r) is a simple way to start the analysis. 
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5.1 Possible generalizations to 0(r) and 0(s,r). 

Static spherically symmetric solutions. 

Consider first spherically symmetric solutions depending on all 3 coordinates, i.e. 
on r = + z'^. Moreover, generalize to n space dimensions. The Lagrangean 

(2.1) becomes 

c = rV-i [l - v/l + /^(0'2 + m202)l ^ (5.1) 

and its equation of motion is 

+ - —^0' — [1 + + m^0) = 0. (5.2) 


where ' denotes differentiation with respect to r. After simplihcations, it is rewritten 
as 

n — 1 _ f „ n — 1 \ 

(5.3) 


ij' -I- - 1 ( 5 ' _ ,7^2 < ^ ^4 


r 1 + 

One dimensional solution. 

In n = 1 space dimension, the ansatz 

Hr) = 

1 + /3r 


2m^0 — -—1^' 


(5.4) 


is an approximate solution. Indeed upon substituting this ansatz into the equation 
of motion, we obtain 


202 -m2(l + 0r)2 = 


202 


1 + 


(l+/?r)^ 


(5.5) 


after simplifying by a common factor A/(l + 0r)^. The ansatz satishes the equation 
of motion, if 0 3> m and 

202 


= m , 


(5.6) 


and then it is valid even in the 0r ~ C^(l) regime, since then we can approximate 




. (l+/3r)2 




(5.7) 


In conclusion, the solution is 


0(r) ~ 


A 


1 + 


Ar 1 
\/2m'^P 


(5.8) 


and as we can see, it is parametrized by A, and is valid for A 3> . 
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However, this solution is only valid in n = 1 space dimension. 

At very large distances, the wavefunction in any dimension becomes 


0(r) = He-™" , (5.9) 

which is what Heisenberg considered as well. 

Solutions of the form (j){s,r). 

A possible generalization that would include both the 0(s) near s = 0 and the 
0(r) near r —)■ cx) is 0(s,r). Substituting this ansatz in the DBI action, we obtain 

2 


hrst 


= -4s ( ^ ) + 


ds 


k/2 


(5.10) 


and then for the Lagrangean 


C = /-V 




\ 


1 + \ — As [ -f- ] + 

ds 


(5.11) 


where n is now the number of transverse space dimensions (n = 2 in the physical 
case). As before, we hnd the equation of motion 


t ft ^ ^ / 9 / A ^ 

h H-0 — — 4— 

ds 


r 

+8Hs( ^ 

ds 


dcj) 


dcj) 

— + m 
ds 


n — 1 


•A A/ 


dcj) 


ds 


ds 

2 “ 


[1 + 


4s/N^) 0"-4/V'^ 


ds 


ds 


d(j) 

ds 


4 d<t) ^,d<t)' 


+ = 0 , 
ds ds 


(5.12) 


where the third line contains terms with mixed derivatives. 


We can check that again at s ~ 0, 0 = A^/s is a solution, but 0 = Ay/sf{r) is 
not a solution at nonzero r, since the leading terms in the equations of motion for 
such an ansatz are 


0 ~ —4-^ 
ds 


dcj) 

ds 


d(l)\ 


A 


A 


+ 8i<q^l =-^f(r) + ^AVn{r), 


(5.13) 


/—j j ' r 

V-s 

and this equation has as only solution /(r) = ±1. Therefore the solution at nonzero 
r and s ~ 0 must be of the type 


0 AaTs + s'^f{r) , 


(5.14) 


where n > 1. 
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5.2 Delta function shockwave? 


Before we continue with 0(s, r), we want to address the issue of a possible delta func¬ 
tion shockwave. In the gravity dual theory, the gravitational shockwaves that scatter 
are delta function shockwaves [9, 10], so one can ask whether the same happens also 
in the held theory picture. 

We want then to try a delta function ansatz for a 0 = (j){x~,r), where x~ = 

{x - t)/V2, 

(j){x~,r) = 6{x~)^{r). (5.15) 


The equations of motion for (j){x , r) in n = 2 transverse dimensions are 



r 


k/2 




1 -I- 




(5.16) 


Then for the delta function ansatz, in the denominator on the right hand side of 
(5.16) we have the 1 negligible with respect to the 0^ term (since the whole term 
is proportional to a delta function, so the denominator is relevant only on the delta 
function, when the value is inhnite), and the equation becomes (after simplifying the 
common delta function on both sides) 


-h -<h' - 
r 


^$'2 / ^ 

$ \ 2m‘^r <l> / 


(5.17) 


But we can easily verify that this equation has no solutions of the form at 

large distances, nor of Ar~P type, and if we put Ae~°'^^ we hud that the only solution 
is p = 2 , 0 ; = —m^/2, i.e. which is clearly nonphysical. 

We do have in fact the solution Ae''^^, but it is a complex solution for a real 
scalar, and Acos{mr) is not a solution (the equation is nonlinear, so we do not have a 
superposition principle). So the conclusion seems to be that this case 0 = (5(a;“)<h(r) 
is unphysical. 

In fact, there are some ways around that. We can consider a case when the 0 
resembles much a delta function, but it has a hnite thickness, and the height of the 
delta function is not only hnite, but such that the 1 in the denominator of (5.16) 
actually dominates, so we get the equation of motion 

(^2m$ - . (5.18) 


Another possibility is to add by hand a source (pS^x )/(r) to the action, leading 
to the modihed equation of motion 


1 *'2 

= 2 — ( 1 - 

r $ 


1 

2m^r <l> 


$'2 \ T 2 


(5.19) 
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but we are still left with the issue of understanding the source-free shockwaves like 
Heisenberg’s. 

Instead, we can notice that we do not really need a delta function shockwave in 
, only need that T+_|_ becomes inhnite at a;’*' = 0. Normally that happens because 
of a in T+_|_ which implies also a in the field (in the case of the gravity 

dual, delta function in the metric). But in the case of the solution of Heisenberg, we 
just have an energy density that blows up slowly near x'^ = 0. 

Indeed, near s = 0, we have (see (2.11)) 






ruy/s 


I '^x'^ 
3J2 


ms 


(5.20) 


which blows up at s = 0. Moreover, we can calculate T++, which turns the in the 
numerator into s, implying T+_|_ ~ l~'^/my/s —)■ 0. That means that there is a source 
at s = 0, since T+_|_ becomes inhnite there, i.e. at a;’*' = 0 and = 0 (two plane 
waves, travelling in opposite directions). In the next section we will study the source 
of the pion held in more detail. 


5.3 Perturbative solution near r = 0 


We now return to 4>{s,r) and consider the expansion near r = 0 of (j){s,r). We have 
fonnd the equation of motion (5.12), and the ansatz (5.14). We hrst plug this ansatz 
in the equation of motion for n = 1, bnt we hnd that while at zeroth order we get 
zero for A = then at hrst order we do not have cancellation. 

It means that we need to consider the next order in ^/s, namely n = 3/2. Then 
we can check that the relevant terms are only 


+8l‘^s I 


d' 

ds J 







-|- QA?l^f \/~s -|-... , 


(5.21) 


but now we see that with A = we cancel both zeroth order and hrst order 
terms. Moreover, now we have other terms in the equation of motion contributing, 
in particular 


2 / A ^ 

—m(j) — 4— 

ds 

= —m^A^/s — A 


' d 


I m (f) + 81 sm (j) 


d(f)\ 
ds J 
0 , 


2 


(5.22) 


but these also cancel! 
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That means that we need to consider also the second subleading term in the 
expansion in s of (j){s,r), 

0 = A\/s + + s^^‘^g{r) , (5.23) 

and check the terms of order in the equation of motion as well. Incidentally, we 
can check that considering a power s°‘ smaller than 5/2 doesn’t work either, since 
then again only the two terms above contribute to second subleading order, the hrst 
giving —4a^g{r)s°‘~^, and the second giving +Qag{r)s°‘~^, so they only cancel for 
a = 3/2, which is excluded (is the hrst subleading term). 

Then we obtain 


d 

ds [_ ds 



-9fVs -25gs^/^ 


V~s 


+ 9A^l*f^ + 27 + 15kl2/V^/2 + ... (5.24) 


The other terms on the hrst line of (5.12) give to order 

^ 3/2 j-rr ^ _ ^2 j _ ^ (5.25) 

r J 

the terms on the second line do not contribute to this order, and the terms on the 
third line give Summing up all the contributions, and using the zeroth 

order condition A = we obtain 

(5.26) 


^3/2 27/^/^ — 19g + - — -f' + AQm^f — m^l 

r 


and equating this to zero hxes gigr) to be 

1 79 — 1 

gir) = — 27/^/^(r) H-/'(r) + 46m^/(r) — . (5.27) 

10 L r J 


The interpretation is that we can specify arbitrarily the function /(r), or in 
another way specify the function 



(5.28) 


which is an initial data on the Cauchy surface s = 0. Once this is given, the rest of 
the function 0 should be hxed by the equation of motion. 
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5.4 Perturbative solution near r = oo 


However, we are interested instead on the behaviour at large, but hnite r, needed for 
the calculation of cTtot('S) through femax- 

We know that at small held and derivatives, the DBI action reduces to the free 
massive scalar action. Indeed, viewed as an expansion in /^, or in nonlinearities of 
the held, the equations of motion reduce to zeroth order to the free equation 


^ ~ 2 , . ^ 

H-0 — m 0 —4— 

r as 


dcj) 

ds 


= 0 , 


(5.29) 


and so, under the assumption that the s dependence is subleading, and we can ignore 
the last term involving only d/ds, we obtain at mr 3> 1 the solution 


0 




Ae 


(5.30) 


as expected. Note that the 0' term in the equations of motion is subleading in mr 
and it doesn’t contribute to this order. 

But we can be more precise, since the exact solution to the free equation of motion 
is known. If we had n = 3, the exact solution would be the Yukawa potential, 

Ap-mr 

0(r) = 5Y-. (5.31) 

r 


For n = 2, the exact solution is a bit more complicated. We can write the 
equation of motion at nonzero mass m as 


d^(j) 

d{imry 


+ 


1 d 
imr d{imr) 


0 + 0 


0 , 


(5.32) 


which matches the dehning diherential equations of the Bessel functions at index 

z/ = 0, 


dz"^ 


+ 


1 dZ^ 
z dz 


+ Zq — 0 , 


(5.33) 


and therefore we have 


0 = Z^iimr). 


(5.34) 


We want to choose the Bessel function of imaginary argument that decays expo¬ 
nentially at inhnity. This is 

K^imr) = , (5.35) 

giving for the scalar 

0(r) = AKo{mr). (5.36) 
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The asymptotics at mr —)■ cx) give 


(f){r) ~ A 


TT 

2mr 


but we should also note the asymptotics at mr —)■ 0, where 


Z Z 

Kq{z) ~ - In ~ - In 


(5.37) 


(5.38) 


To hnd corrections to this free solution, we could think of expanding the equations 
of motion in or equivalently the mass dimension of the remaining expression (once 
is removed), but besides the free terms above, all the other terms are linear in 
that is, of mass dimension 7 with respect to the rest. 

We can instead take an ansatz that 0 depends only on r, and not on s, with a 
coefficient that is of the order of i.e. 


(f) = AK^imr) + Bg{r) , 


(5.39) 


and B Then the full equation of motion reduces to 


B ( g''{r) + -g'{r) - m^g{r) ) = 


- (f)' 


(5.40) 


With the assumption that B oc we can consider on the right hand side only 
the order zero term with A, to obtain 


B (^g'\r) + ^g\r) - m^g{r)^ 
= {mlYA^ ( ^ Ko{mr 


d{mr) 


2KQ{mr) 


d 


mr d{mr 


-K^imr) 


(5.41) 


However, even for that, we can only hnd the leading order solution at mr —)■ cx). 
Then g{r) ~ e“^™^/(mr)^/^ solves the equation to leading order, and we hnd the 
solution 

'7r\3/2 g-3mr 


0 ~ AKoimr) + j 


{mr^B 


(5.42) 


at mr —)■ cx. 


6 The source for the pion field 

In this section we would like to understand what is the source of the pion held, 
which is supposed to represent the nucleons. To do so, we hrst look at the original 
Born-Infeld action for nonlinear electrodynamics, in terms of a held strength 
and apply the lessons learned to our DBI case, hrst for a static solution, then for the 
shockwave. 



6.1 The vector Born-Infeld case 


In the original paper of Born and Infeld on nonlinear electrodynamics [19], the issue 
of the source for solutions of the nonlinear Maxwell held was explored, and in fact it 
was the crucial motivation for the work: to obtain a smooth ’’electron” solution to 
the equations of motion, free of singularities. 

The BI Lagrangean can be written as 

c = VTTy^^-I , ( 6 . 1 ) 


where we dehned 

G^^iB-E). (6.2) 

Using these dehnitions we can dehne quantities analogous to the quantities dehned 
for electromagnetism in a medium, namely 


H 

= = 

B-GE 

dS 

Vl + E-G"^ 

D 

= = 

E-GB 

dE 

Vl + E-G^ 


(6.3) 


The equations of motion and Bianchi identities of the BI Lagrangean, that correspond 
to Maxwell’s equations of the linear theory, are 


V X E + doB = 0; V ■ .B = 0 

W xH-doD = 0; V ■ D = 0. (6.4) 


The Hamiltonian density can be written as 

H = Vi + p - _ 1, (6.5) 

where 

^ = Q = y2^d.h), 

The inverse relations for the helds are obtained from the Hamiltonian as 

5 H + QD 

dH v^i + p - g2 

^ d + qh 

dD v^i + p - g2 

At zero magnetic held, B = H = 0, the equations of motion and Bianchi identities 
reduce to 

VxP = 0; V-P = 0, (6.8) 


( 6 . 6 ) 


(6.7) 
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and we also find Q = 0, P = D'^/b'^, G = 0, F = —E‘^/b‘^. Then V ■ -D = 0 reduces to 

^y-D,) = 0, ( 6 . 9 ) 

which admits a non-trivial solution of the form 

D. = ^ , (6,10) 

the same as in Maxwell theory. More precisely, the solution is at r 7^ 0, which means 
that we have actually 

V ■ D = A'Ke5'^{r) , (6-11) 

which gives the integral formula (from Gauss’s law) 


47re 



( 6 . 12 ) 


So from the point of view of D (the field in the medium in electromagnetism, 
where D = e^E + P and V ■ T* = 0 in the absence of external sources, and otherwise 
just includes the charges external to the medium) the sources are point-like as in 
Maxwell theory. 

For a static system we have E = — VAq, where Aq is the zero’s component of the 
gauge held vector potential, and it is related to D by 


D = 


E 


— = Dr- = 


Er- 




b2 


1 — ^ \ 1 — — 


(6.13) 


which implies for the electric held E 


— Er — A!q — 


e/rl 


+ D/rl ’ 


where 


ro = M-^ 


(6.14) 


(6.15) 


is a radius related to the radius of the electron. We then also obtain the electric 
potential 

p / r \ 

(6.16) 


Mr) = —/ — , 
ro V^o/ 


where 


fix) = 


dy 




(6.17) 
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and we obtain that /(O) ~ 1.8541 and Aq ~ 1.8541e/ro. 

The hnite maximum of the electric held E is obtained at r = 0, and equals 
e/rl = b, as expected, since C = ^1 — E‘^/b‘^. 

As in electromagnetism, we can regard the total charge Q in the material as either 
J dW ■ D or J dW ■ E with the difference that the former expression counts only 
the outside charge introduced, whereas the latter expression counts all the charge, 
including the polarization response of the material, which tends to spread out the 
charge density. 

Also now, we can define 47rp = V ■ E and find after an easy calculation that 


e 1 

^ 27rr^ (r/ro)(l + (r/ro)"^)^/^ 


(6.18) 


and we see that this charge density is spread out, going as 1/r^ at r —)■ cx), but only 
as 1/r at r —)■ 0. We can also verify the fact that its integral gives the same result 
as the integral of D^, namely e. 


6.2 The scalar DBI action and its source 


Static scalar DBI results 

A similar thing happens for the scalar DBI action. We start by reviewing the 
construction of the static scalar solutions paralleling the nonlinear electrodynamics 
solutions, as presented in [20]. 

On static solutions, dt(j) = 0, the scalar DBI action reduces to 

£ = ^1 + F2 , (6.19) 

where 

F = V(j), (6.20) 

0 being the DBI scalar. Note then that this action is the same as the vector BI action 
above for the case B = 0, just changing the sign inside the square root. Therefore 
we can follow the same analysis, and first define 


^ dC F 

dF Vi + F2 

In terms of it, the equation of motion is 

V-O = 0 , 


( 6 . 21 ) 


( 6 . 22 ) 
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solved by 




(6.23) 


so that really we have V ■ C* = 47re in 3 spatial dimensions. Therefore the solution 
for the scalar is given by 

Fr = dr(p = , , (6.24) 

v /1 - e^lr^ 

which is called the ’’catenoid”. The solution has a horizon-like structure at r = 
^Je = tq, due to the fact that it has the interpretation (in the case it is the action of 
a D-brane) of one half of a D-brane-anti-D-brane solution connected by a throat. 

The DBI scalar shockwave- 

In the spirit of the model of [3] we now consider a four dimensional scalar held 
0(r, s). In particular (j){s) can be recast from a 1-1-1 dimensional action which for the 
massless case reads 


C = I 


-4 


as 


(6.25) 


We dehne hrst the analog of the electric held from the Born Infeld paper, 

Es = 2^^. (6.26) 

as 


In terms of it, the Lagrangean becomes 

C = l-^ 


(6.27) 


just like the BI vector case. Then we also dehne the analog of the electric induction, 

dC 


D. = 


which gives 


D.. = 


dEs ’ 

E.. 


The equation of motion (the analog of Maxwell’s equation) is 


solved by 




Ds = s> 0. 

Vs 


(6.28) 


(6.29) 


(6.30) 


(6.31) 
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Causality then requires that we have Dg = 0 for s < 0. Then inverting Ds{Es) we 
get 

Eg = , ■ (6.32) 

But since = 0 for s < 0, we also have Eg = {) for s < 0, which means that really, 


P - ^ PIP- F - 
Es Eg . - . „ • 

V-s V-s + 

We can also integrate the above to find that 0 is given by 


(6.33) 


= [ ds^ = ^ [ ds 

J 2^ 2 J 


a / s(s + EA^) 


= A log 


a/s + a/s + EA'^ 


(6.34) 


at s > 0 and 0 at s < 0 which has the same structure as of (2.8). 
This reduces at small s to 


0(s) ~ I ‘^^/s9{s) , 


(6.35) 


which is the same solution as Heisenberg’s. Note that the constant A determining 
Dg is arbitrary, even though 0(s) near s = 0 is completely determined. 

Then the electric field is a step function. 

Eg = 2^/:s^ ~ l-‘^eis) , (6.36) 

ds 

and the electric induction is 

Dg = (6.37) 


Plugging back in the equation of motion for Dg, we have really for the analog of 

Pext. = V ■ D, 

^iV^sDg) = UAe{s)) = A5{s). (6.38) 

ds ds 

So as in the BI case, there is a source term, which is a delta function when 
viewed from the point of view of the induction Dg (i.e., it is an ’’external source” to 
the medium). The value of the charge. A, is arbitrary, even though 0(s) near s = 0 
is completely determined. 

We can also dehne the equivalent of the V ■ E = p, the total charge (including 
the one due to the ’’polarization of the medium”), which is spread out. We dehne 
the density 


P 


A 

ds 




( V~s0{s) 5 ^ 

\^s + DAy 


2v^(s + H/M2)3/2 


(6.39) 
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Note that we dropped a term coming from the derivative of 0{s), proportional to 
■\/s(5(s), since this is zero. We see that this charge drops at inhnity as 1/s^, and at 
0 only as l/^/s, and integrates to the same total value as the one dehned via Dg, 


A 


V^0(s) 

Vs + /M2 


CXD 


0 


A. 


(6.40) 


In conclusion, there is an ’’external source” located at s = 0 (the shock’s position), 
with an arbitrary charge, but the ”in medium” source is spread out, over an s of the 
order of VA^. 


7 The cross section 

We can now hnally consider the calculation of cross sections arising from the Heisen¬ 
berg model. 


7.1 Corrections away from the Froissart limit 


The hrst issue to address is of a systematic expansion away from the limit of Froissart 
bound saturation. It is clear that by considering a (j){r) that is not yet completely 
dominated by the e"™” term, we can hnd corrections to the Froissart behaviour of 
the cross section. If we have an exact wavefunction, we can obtain a cTtot('S) that 
would be different in the leading behaviour, like a power law atot(5) oc appearing 
before the onset of Froissart saturation. 


Corrections to leading behaviour 

We hrst consider corrections to the e”™” behaviour of 4>{r), which were found 
in (5.42), with the free part being asymptotically (5.38). The behaviour is 

subleading with respect to the 1 / V^ in the hrst factor, so we consider 


^(r) oc 


'mr 


(7.1) 


Then as usual, the emitted energy is proportional to 0(6) V^, and when it gets down 
to {koV (the average per pion emitted energy), we reach 6max- Thus 


Vb 


max'^TT 


- {ko,V , 


giving 


6max ^ - In 


~ -In 


(E^)^ln(v^/(V)) """ 


In 


In 


(^0,7r) 


(7.2) 


(7.3) 
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and atot(s) = T^b max(^) • 

Possible new regime 

But besides the small corrections to the Froissart saturation regime above, we 
can in principle have also a situation where a new regime for (Jtot('S) appears. 

To avoid the leading Froissart behaviour, we must avoid the exponential for 
f' = &max, so we need to have m^6max('S) < 1. This can indeed exist in some energy 
regime s, for small mass m = 

Since the scale I in Heisenberg’s DBI action can presumably be identified with 
^QCD, and Aqcd ~ 2m^, the corrections of order (m/)^ ~ (^^tt/Aqcd)^ = 1/4 are 
small, so it could be a good approximation. 

But if ml 1, there is a regime where the wavefunction is linear, and when 
solving for (j){r) from the equation of motion we never get into the nonlinear regime. 
That means that the full solution to the free equation, 0 = AKo{mr), is exact. At 
distances r m~^, we obtain 

HIT 

0(r)~—Aln-^. (7.4) 

Then the condition for b^ax at energies s for which the above 0(r) are still in the 
linear regime is 

giving 

1 _2Ashlll 

bmaxis) = -e ^ ^ 

n 71 —aAOiIlI 

O'tot(5) = 7rbraax{S) = , (7.6) 

for > (/cq.tt), which gives a mildly increasing dependence, that could be easily 
mistaken for a small power law or the log^ behaviour of Froissart saturation. 

In conclusion, such a new energy regime could appear in QCD just before the 
onset of Froissart saturation, but it would be hard to distinguish experimentally from 
the small power law (’’soft Pomeron”) behaviour, or from the Froissart saturation 
behaviour. 

7.2 Black disk model and ratio of elastic to total cross sec¬ 
tion 

Until now we have discussed the total cross section, or in the case of several mesons, 
also individual meson cross sections. But we want now to discuss also the elastic 


'2 l^(^'7r^max('S)) 




(7.5) 
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cross section. For that however, we need a quantum amplitude, whose forward part 
gives the total cross section, and whose absolute value squared gives the elastic cross 
section. 


Since we do not have a quantum amplitude, only a total cross section, we can 
engineer an amplitude that gives this total cross section, and from it calculate the 
elastic amplitude. The simplest model is a black disk eikonal amplitude, with S 
matrix S = and Im((5) = cx) for 6 < b^ax{s) and with <5 = 0 for 6 > 6max(5). This 
reproduces the cross section 7r6max(-s)^. 

For massless states scattering, we have in general 


-^(s, t) 
s 


-i j - l) 

/*^max(^) /‘27r 


bdb 




COS 6 




(7.7) 


where b is the impact parameter (transversal), and its Fourier conjugate is q, with 
^ = t. 

For the black disk eikonal. 


1 bi'<- 

S y/t 


-Ji{Vib 

max (3)). 


(7.8) 


The total cross section is found from 

^ImAiastic(s, t = 0) = atot{.ki, -t anything) , (7.9) 

s 

and it is easy to calculate that for the black disk eikonal we get atot = 7r&max('S)^. 

We should note here that most of the times, like for instance in [21], the black disk 
eikonal model starts with a partial amplitude ai{k) = — l)/(2i), suggested 

by the partial wave expansion, which is a factor of 2 smaller than (7.7). After 
the normalization of the cross section is properly taken into account, this leads to 
(Ttot = since atot ~ Ima, but del ~ so a rescaling of a leads to 

a rescaling of dei/dtot, to an aei = TTfe^ax- But our model, also used for instance 
in [9], is physically different, since we considered simply, as usual, the amplitude 
as the Fourier transform of the T-matrix, and A = 1 + iT = e*^. This leads to 
o'tot = 7r^max('5)) which we believe is a model more deserving of the name black disk, 
as the total cross section equals the classical one. Then, as we shall see, we obtain 
o'ei/o'tot — 1/4, instead of 1/2. 

In the case that the particles are massive with mass m instead, the 1/s is replaced 
by 1/ {2pcmEcm)- But if mi = m2 = m, Ecm = 2 a/+ m^ = y/S, so we have 

2 EcmPcm = a/s(s - 4m2). (7.10) 
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It remains to define the range of integration for t, given s. In the center of mass 
system, 

■5 = ^ = {PCM — kcM^ = + PcM ~ ‘^kcMPCM COS 9 , (7-15) 

where pcM and kcM are momenta of the same particle, before and after the collision 
in the center of mass. Then the range of integration for t, given s, which fixes pcM 
and kcMi is 

f ^ [{PcM — kcuY 1 {pcM + (7-16) 

But pcM = kcM and V^/2 = Eqm — E — a/ Pcm + m^, meaning that 

Pcm = kcM = ’ (7-17) 

and then the range of integration of t is 

te[0,s-4m2], (7.18) 

so that hnally 

-1 PS—Am? Jj. 

= i / y[Ji(v^t„„(S))|^ (7.19) 

o-tot 4 Jo t 

By using the recurrence relations for Ju{x), we do the integral and obtain 

^dastic ^ 1 1- (^Jo(|6max(5)|Vs-4m2)j - Jl (| &max(5) | Vs - 4m2) j . (7.20) 
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At large 



1 -^ 

^^max(5) 

where we put rriN for a nucleon or nucleus mass, corresponding to the case when we 
collide nucleons or nuclei. Then from the Heisenberg model &max(>s) — 1/m.jr ln(s/so), 
so that the sought-for ratio is 



^elastic ^ 
_ _ 

<^tot ~ 4 


^elastic 1 ^Tfl^ 

- ~ — 1 - - 

o'tot 4 7rln(s/so)\/5 — 4m^ 


(7.23) 


asymptoting very fast to 1/4. 

This compares very well with the experimental results from the TOTEM experi¬ 
ment [22]. 


8 Heisenberg model and Holography 

In section 3.3 we described a sigma model in AdS space. This can be directly related 
to another holographic model, the “hard wall” model, which is an AdS background 
chopped off at a certain value of the radial coordinate. This scenario is addresssed 
in the following subsection. We then present an alternative approach that includes a 
systematical analysis of the relations between Heisenberg’s model and the holographic 
description of a proton-proton scattering. This in fact involves two steps. In the first 
we will establish the relations between the DBI action used in Heisenberg’s model and 
the DBI action that emerges as the action of flavor branes in conhning backgrounds. 
The second step is to layout the holographic dual of scattering of baryons and to 
relate it to the extraction of the cross section from Heisenberg’s model. The two 
steps are described in the second and third subsections of this section. 

8.1 The relation to the holographic “hard wall” model 

The remarkable fact is that, even though the Heisenberg’s model was proposed before 
string theory was discovered, the DBI action used by Heisenberg emerges naturally 
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in holographic models of QCD since it relates to the effective action of open strings. 
In the simplest model for high energy QCD scattering introduced by Polchinski and 
Strassler, one considers an AdS^ space, 

dT‘^ 

ds^ =—dx^ + , ( 8 . 1 ) 

cut off at an rmin = R^A, with A identihed with the (pure) QCD scale (glueball scale). 
It was soon realized that one can think of the IR cut-off as a dynamical IR brane 
(like in the Randall-Sundrum model), and the appropriately normalized fluctuation 
in the position rmin, the scalar ~ 0 can be identihed with the pion in QCD. But the 
action for the huctuation in position of a brane is exactly the DBI action! 

The only nontrivial part of the action is the potential for the brane position, 
which can appear, depending on the mechanism, either inside or outside the square 
root. 

The picture of high energy scattering is also similar in the gravity dual [9, 10, 
11, 12, 13, 14]. In a purely gravitational theory, we have gravitational shockwave 
collision, happening near the IR cut-off, creating a black hole on the IR cut-off, being 
mapped to the pion held shockwave collisions creating a hreball. But more precisely, 
when we consider also the huctuation of the IR cut-oh giving the pion, we have the 
same picture, of pion held shockwaves colliding and creating a hreball. 

8.2 The DBI action of flavor branes in conflning backgrounds 

Heisenberg’s model assumes that the scalar helds that are in charge of the interac¬ 
tion between nucleons are governed by a DBI action in hat space-time. Holography 
provides dual string descriptions to certain strongly coupled gauge dynamical sys¬ 
tems. As was mentioned above, the DBI action is a basic tool in the toolkit of string 
theories. Thus, an obvious question to ask is whether one can relate Heisenberg’s 
model to a holographic description of proton-proton scattering, and in what way. 
To answer this question one has to address hrst the issue of what is the holographic 
laboratory dual of QCD in its conhning phase. 

The basic AdS^ x string theory, the dual of A/” = 4 SYM is clearly not the right 
setup. It is both conformal and maximally supersymmetric. One has to deform the 
geometrical background in such a way that the isometry group is not SO {4, 2) x SO {6) 
but rather only the four dimensional Poincare symmetry. Obviously the desired 
background should be equipped with a scale which breaks scale invariance. To check 
whether a given background corresponds to a boundary conhning held theory, one 
should investigate the stringy dual of the Wilson line. A necessary condition for 
a “conhning background” is that any rectangular Wilson line along a space and 
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the time directions should admit a conhning area law behavior. In ref. [23] it 
was shown that this is achieved provided that either gttdxxiu) has a non-vanishing 
minimum value or that it does not vanish at the value of the radial coordinate u 
where gudnuiu) —)■ oo, and where gtt,gxx and guu are the metric components along 
time, the space direction of the Wilson line, and the radial direction respectively. 
Not surprisingly the AdS^ x background does not obey this requirement. 

A close cousin of this background that does admit conhnement is the “hard wall 
model” discussed in the previous subsection where one, by hand, chops off the radial 
direction to be m > u\ where ua is a scale in the bulk that corresponds to Aqcd of 
the boundary conhning gauge theory. This however is not a solution of the equations 
of motion. 

A prototype conhning background that is a solution is the AdS^ background 
with one spatial coordinate compactihed [24] on a circle in such a way that the two 
dimensional manifold spanned by the radial direction and the circle has a cigar-like 
geometry. It is easy to check that upon imposing anti-periodic boundary conditions 
for fermions, the only massless helds of the dual large W gauge theory are only 
the gauge helds and all their supersymmetric partners become massive. In that 
way supersymmetry is broken and the dual held theory is that of pure large W 
gauge theory in three space-time dimensions. To get a gravity model dual of four 
dimensional conhning large W gauge theory, one can compactify the near horizon 
background of large number of DA branes [25] rather than the AdSjj x model which 
is the background of large number of D3 branes. In fact the dual gauge theory is 
an ehective conhning theory with energies smaller than ^ where R is the radius of 
the compact circle which maps into the mass of the dual glueballs. There are several 
other solutions of the ten dimensional supergravity equations of motion that admit 
conhnement but with no loss of generality we will discuss here only this model. 

To incorporate in the gravity side the quark degrees of freedom one introduces Nf 
havor D-branes. For Nf << Nc one can neglect the back-reaction of the havor brane 
on the bulk and hence treat them as probe branes. In the Sakai Sugimoto model [15] 
a stack on Nf D 8 and a stack of anti- D 8 branes are placed so that asymptotically 
at large radial direction the transverse direction to their worldvolumes is along the 
compact circle X 4 . 

In the IR in the region of the tip of the cigar the two stacks of branes have to merge 
one into the other hence breaking the original UiiNf) x Uii{Nf) chiral symmetry into 
a diagonal subgroup of UoiNf). In original model, the U-shaped branes were in an 
antipodal setup uo = ma , see the right figure of ( 2 ). This was generalized (see the 
left figure) to incorporate an additional parameter uq 7 ^ ua [16] which, as will be 
shown below,is crucial for coupling the protons to pions in the holographic picture. 
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Figure 2: On the right side we have the antipodal set-up of the Sakai-Sugimoto model 
where Uq = u\. On the left side we have the generalized non-antipodal set-up. 

The physics of the degrees of freedom that resides on the flavor branes, namely the 
U{Nf) gauge helds and the scalars in the adjoint of the U{Nf) group, is described 
by a DBI action. In fact the action includes, on top of the DBI action, also a CS 
term. That is obviously where Heisenberg’s model and holography meet. The action 
on the flavor branes in the Sakai-Sugimoto model reads 

Sdbi = Ts f dVe-^- det[d^X^d,X^gij{X) + 2TTa'F^,] , (8.2) 
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where the dilaton 0, the metric Qij and the RR four form are given by [16] 


ds'^ = 


—) 

Rd4 ) 


3/2 


-dt‘^ + 5ijdx^dx^ + f{u)dx 


+ 


(^) 


3/2 


du^ 


Ifiu) 


+u^dQ.l 


(8.3) 


F 27riVc ^ 

Fa = -^£4 , = gs 

Va 


— ) 

Rda ) 


3/4 


. = ^9,NX. , /(«) = 1 - ( — 

u 


where xa is the coordinate of the compactihed circle, Va is the volume of the unit 
four sphere 114 and 64 its corresponding volume form. Upon inserting the metric and 
the dilaton one hnds, according to the general analysis in section 4.2, 


Sdbi = Tg / dtd^xdxA 0^i//(0) + 


R 


DA 


3 r 




/( 0 ) 


{dxM' 


(8.4) 


where Tg = TgQ^/gs and to connect to the rest of the paper we denoted the radial 
coordinate u by 0. 

The fluctuations of 0 translate using the dictionary of holography to scalar 
mesons. To extract the spectrum of the latter one considers hrst a prohle of the 
flavor brane given by (pciixA)- One then introduces the fluctuations of 0 in the fol¬ 
lowing form 

0(0:4, X^) = 4>ci{Xa) + d4>n{x^)Cn{xA)- (8.5) 

n 

The lowest mode of the fluctuating held 0o should be identihed with the scalar held 
(j){x^) in Heisenberg model. Next one expands the DBI action to quadratic order in 
0, integrates over the 0:4 direction, and derives a massive spectrum for the (50„(a:^). 
Here we do not want to expand the square root but rather maintain the full tower 
of derivatives of the held. The outcome of the integration of the Cn(o:4) will be 
mass terms of the form m^0^ plus terms higher order in 0. We assume here that 
the truncation to only the mass term in the expansion of the DBI action can be 
translated to having a mass term in the four dimenional DBI itself. In that case the 
action takes the form 


Sdbi — Tg 


dtd^x(l)'^\ /(0) -|- 



[(9^00^0 -f 771^0^ . 


( 8 . 6 ) 


The equation of motion that associates with the action (8.6) for the m = 0 case can 
be written in the following form 




8-5 


+8 


Rda 


.u\R%a 


3 / D \ ^ 

Kda 


(( 9 ^ 0 ( 9'"0 - 200 ^( 9 ^ 0 ) 


\2l2 


d,<pdy - - 2 ^[( 0 , 0 )^] 
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+2 


6 


(8.7) 


R 


Di 




We leave the investigation of the relation between the solution of the DBI action 
given here and the DBI in flat spacetime used in Heisenberg model to a future 
investigation. 


8.3 A holographic description of the proton-proton scatter¬ 
ing 

So far we have discussed the holographic laboratory and its relation to the DBI 
action applied in Heisenberg’s model. Next we would like to see what is the relation 
between the cross section of a proton-proton scattering in Heisenberg’s model and 
the corresponding cross section in a holographic setup that associate with a conhning 
theory equipped with flavor degress of freedom. Here for concreteness we will use the 
generalized Sakai-Sugimoto model. A stringy realization of a baryon in this model 
[26] is that of a baryonic vertex made out of a D4 brane that wraps the four cycle and 
is connected by W strings to the Nf probe flavor branes [27]. A priori the baryonic 
vertex could have been located in the generalized Sakai Sugimoto model in any place 
below the flavor brane, but in [28] it was shown that in fact it must be immersed on 
the flavor brane. The interaction between two protons in this setup is that of two 
baryonic vertices each connected to Nc strings that stretch on the flavor branes. The 
scattering of such two objects is obviously very complicated. Instead it was shown 
in [17] that one can view the baryon as a flavored gauge instanton. This follows 
from the fact that the wrapped D4 brane is a point on the four dimensional part of 
the world-volume of the flavor brane which is spanned by the ordinary three space 
coordinates and the radial direction. Alternatively it can be shown by expanding 
the flavor gauge DBI-I- CS actions, keeping only the leading order U{Nf) YM -|- CS 
action. The 5 dimensional action takes the form 


S = Sym + Scs 


Sym ~ 
Scs = 


' crx / dz 






Nc f 2 ^2 f . ^2 

‘ AA trF^ + ^ I A AF\ 


where near the bottom of the U-shaped flavor branes we have 


1 


_ _ NcXMkk 

~ 2167r3 


C + + 0{A4^z^: 


( 8 . 8 ) 


(8.9) 
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Here T is the U{Nf) gauge field, and A and A denote the gauge one-form as¬ 
sociated with the SU{Nf) and U{1) subgroups respectively. We made a coordi¬ 
nate transformation from {x^,X 4 ) to to a five-dimensional conformal metric ds"^ = 



Based on this action it was further shown 

UA 


that the static properties extracted from this model are similar to those derived from 
the Skyrme model. 


Next we would like to examine to what extent does Heisenberg’s treatment of the 
scattering of a proton on proton translate into a scattering process of two instantons 
in the holographic laboratory. The interaction of the latter can be divided into three 
zones[29]. In the far zone when the distances between the two instantons is much 
larger than the inverse of the dual of Aqqd the interaction is dominated by the 
exchange of the lightest meson. In the isoscalar channel it was found out that the 
repulsion, due to the exchange of vector mesons, is stronger than the attraction, 
due to the exchange of scalar mesons, since the lightest meson on the latter type is 
heavier than the lightest vector meson. In the isovector channel it is obvious that 
the lightest meson is the pion and the exchange of it yields an attraction. In the 
near zone, using the solution that carries instanton number equal two, one finds that 
there is only a repulsive hard core interaction. In the intermediate zone there is a 
repulsion due to the interaction of the instanton density with the U{1) of the U{Nf) 
flavor gauge group. However, as was shown in [30] there is also an attractive force 
due to the interaction of the instanton density with the scalar field associated with 
the fluctuation of the D8 branes. The action of this scalar takes the form 


S^ = S 


DBI 


+ 


Nc 

IOtt^ 


U^x dz C{z) X + 


( 8 . 10 ) 


where C{z) measures the ratio of the attractive to the repulsive forces. 


K 

Fr 


= C\z) 


1 - c Y % \ 

9 Yw/ 


< 


1 

9 


< 1 . 


( 8 . 11 ) 


Note that for self-dual (instanton) configurations, Tr[FMNF^^] = Tr[F /\ F] and 
hence the scalar field that originates from the brane fluctuations couples to the 
instanton density, namely to the proton density. 

Thus, in a holographic description of the interaction between two protons, both 
in the intermediate as well as in the far zone, the interaction is mediated by a scalar 
field that is governed by a DBI action. The DBI action (8.4) is not the one Heisenberg 
used but rather a DBI of a scalar in a curved background. The source of the scalar 
field and its coupling to the proton given in (8.10) is different from the source of 
the scalar field discussed in section 5, but a fixed gauge field profile will generate a 
function /(r) in the action as in (5.19), or an implicit external source as in (6.37). 
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9 Summary and open questions 


As was explained in the introduction, in this paper we addressed four aspects of 
Heisenberg’s model of scattering of nucleons: 

(i) We elaborated on, and gave further supporting evidence for the model. We 
made an analysis of the energy of the scalar held, and the conditions under which 
we obtain the (almost) saturation of the Froissart bound. We have analyzed the 
what happens when we go from a 1+1 dimensional solution to a 3+1 dimensional 
one 0(s,r). We have understood the implicit source in the Heisenberg solution by 
analogy with the electromagnetic Born-Infeld action: there is an ’’external” (5(s) 
source that is ’’spread out” by the medium. One can also consider S{x~) shockwave 
solutions by adding an explicit source in the Lagrangean. By using a perturbative 
(j){s,r), we have obtained corrections away from the maximal Froissart saturation 
behaviour, as well as a new regime for atot('S). 

(ii) We examined the uniqueness of the DBI action in terms of giving the (almost) 
saturation of the bound. We have found that, perhaps surprisingly, no action with a 
potential interaction, or with a hnite number of higher derivative terms can do the 
job. The DBI action can do the job, though we have not been able to prove that 
another action with inhnite number of higher derivative terms cannot do as well. 

(hi) Generalizations of the model. We proposed and analyzed several generaliza¬ 
tion of the Heisenberg model. We added a general potential inside the square root, 
instead of just the mass term and we considered a sigma model with several scalar 
mesons. We considered a ’’curved space” generalization inspired by holography, in 
particular the “highly effective action” of [5] for the case of single scalar in AdS^, 
when we replace d^(j)d^(l) by and when considering the power rather 

than (()* we have shown that only for the range n G (0,2) can we obtain saturation 
of the bound. By considering a ’’black disk” type of amplitude in the sense of [9], 
we have obtained also a value for the ratio of the elastic to the total cross section, 
o'ei/o'tot that asymptotically goes to 1/4. We note that the more common model in 
for instance [21] would give 1/2, but the experimental evidence points towards 1/4. 

(iv) We have considered the relation of the Heisenberg model and the DBI action 
he considered to two holographic approaches to proton-proton (or nucleon-nucleon) 
scattering: a simple hard-wall model, and a more precise model based on flavor 
branes in conhning backgrounds. 

In this paper we have just explored the tip of the iceberg. There are a handful of 
additional open questions that are awaiting further investigation. Here we list few 
of them. 
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• Probably the most interesting topic related to realistic high energy scattering 
is performing a precise comparison between the results of Heisenberg’s model 
and experimental data of high energy scattering of nucleons and of nuclei. 
One can deduce the scattering total cross section and the ratio between the 
elastic and total cross sections not only for the asymptotic range of energies 
as was discussed in subsection 7.1. In section 4 we analyzed several possible 
generalizations of the model, and in section 8 we discussed the relation to 
certain holographic models. These deviations from the original model can also 
be confronted with experimental data. One would like to extract the values of 
the parameters of the various models that admit the best £t to the data. In 
particular the mass of the scalar particle that mediates the interaction which 
we referred to as the “pion” in this paper. 

• It is well known that there are two approaches of phenomenological htting the 

experimental data. One is based on the Froissart bound, namely atot log^(s) 
and the other on an exchange of Reggeons and Pomerons between the two 
scattering nucleons. In this case one uses relation of the form atot ~ + 

^g0.08_ Both approaches yield a reasonable £t (see [31] for a possible way to 
connect the gravity dual picture of gravitational shockwave scattering to the 
soft Pomeron behaviour). Thus, a natural question to ask is what is the relation 
between the two models. In section 8 we have attempted to relate the model to 
a holographic model of scattering of nucleons. The latter is an approximated 
picture of a fully stringy description of the scattering process. The exchange 
of a Reggeon and a Pomeron seem closely related to an exchange of an open 
and a closed string. Hence one may be able to hnd a direct relation between 
the two approaches. 

• One natural generalization of the model that was not discussed here but in 
fact is quite common in implementing the DBI action in holography is the 
non-abelian DBI model. To incorporate the (flavor) non-abelian nature of the 
pions is the analog of using Nj probe flavor branes rather than a single one in 
holographic models. In both cases the non-abelianization will provide further 
structure. A first try for the nonabelian model was presented in [20]. 

• Describing the scattering of two nucleons as a scattering of two shockwaves is 
clearly only an approximation and one may attempt at introducing correction 
beyond the shockwave limit. Similarly one can introduce corrections to the 
black disk model. 
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A An alternative method of determining the scalar 
field energy 


The Hamiltonian density was given in (2.11). It reads 


n = p^-c 


+ (V0)2 + 
a/ 1 -h H[((9^0)2 -h 


(A.l) 


To determine the Hamiltonian density in momentnm space namely 'H(A;) is a non¬ 
trivial task for the DBI action since we cannot simply, as is done for ordinary free 
held theories, snbstitnte the Fonrier transform of the held into (A.l) since the helds 
appear also in the denominator. In case that npon snbstitnting the classical solntion 
0(s) into the denominator the latter is a constant then one can use the usual method. 
But as was shown in section 2 this is not the case for the DBI action and hence one 
has to adopt another approach. Here we suggest such an alternative. Dehne now the 
Fourier transform of "H as follows 



dk [e 

2 


(A.2) 


and substitute it into the energy, so that 


E 



■ f [e-^’^^Eik, t) + t)] f [e-^P^E{p, t) + EP^E*{p, t)] 

dX I QjK ^ I CLJ) ^ 

[2E{k,t)E*{k,t) + E{k,t)E{-k,t) + E*{k,t)E*{-k,t)] ^ 


For the theory of a free massless scalar in two space-time dimensions 'H{x,t) = 
\[{dx(t)Y + ('9o0)^]- In this case it is easy to see that E{k,t) = \/ka{k) where the 
held (f>djt) has a Fourier transform (/(x, f) = f dk^^^=[a{k)e~^^^ + {k)e~^^^^]. In 

the case of massive free scalar held we get E{k, t) = \/\/k‘^ + w?a{k). In these cases 
the only contributions to (A.3)will be from the E{k, t)E*{k, t) term. For the general 
case one has to hrst determine E{k,t) and then E is given by (A.3). 

Following this approach we now have to hnd the Fourier transform of 
We cannot hnd an exact analytic expression for it neither for the massless case nor 
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for the massive one. From the analysis of the energy as an integral over s one hnds 
that the main contribution to the energy comes from the region of small s. Thus we 
can get an approximation of the dependence of the energy on 7 using the the leading 
order in s expression of a/'H(s) ~ Fourier transform 

reads 


TF 


\/W) 


(-^)5/8(j274r(|) ■ 


(A.4) 


Expanding this expression in | we get 





(A.5) 


Substituting this expression in the energy and taking the integration region to be 
7 m > k > m we hnally get that 

E ~ ^/ym. (A. 6 ) 

The reason that this result does not match the result found in section 2 is that 
we took a crude approximation of a/'H(s). Obviously this approximation can be 
systematically improved by improving the approximation of 


B Scalar solutions in 0 + 1 dimensions 


Here for completeness we write down solutions of the Heisenberg action in 0+1 
dimensions. The equations of motion in this case are 


0 + m^(j) + 


(j) — 


1 — /[( 0 )^ — 


= 0 . 


(B.l) 


For the massless case the equation reduces to 0 = 0 and hence the solution takes the 
form (j) = at + b. For the massive case the solution takes the form 


y{x) 

y{x) 


isn {im\/lci + lx + im\/lci + 




m 


isn [my/lci + lx + mi/Zci + IC 2 ) | j 


V~i. 


m 


(B.2) 


The solution takes the following form for various values of y = 0.1,1,5 
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Figure 3: 0(t) as a function of t for™ = 0.1 



Figure 4: as a function of f for™ = 1 


For the one dimensional case the Hamiltonian density (2.11) is the Hamiltonian 
and hence we can write a hrst order differential equation which is its conservation in 
time instead of the equation of motion. The Hamiltonian for this case reads 

Hl= ^ :[l + /mV)]-l. (B.3) 

Y1 — i[{^Y ~ 
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Figure 5: 0(t) as a function of t fory = 5 


Thus the first order differential equation is 


= 




(iH + iy 
or in an integral form 


+ m 1 — 




dcj) 


{IH + iyJ^ {IH + iy 


= t + c. 


j H(lH+ 2 ) . 2 (-\ 2 \12 


(B.4) 


(B.5) 


C Scalar solutions in 1+1 dimensions: static and 
depending independently on and x~. 


Before discussing a genuine two dimensional case let’s check the equation for a (soli- 
ton) static solution. For that case the equation takes the form 


dlcf) — — l‘^{dx(j)f 


dl4) + mP'cj) 


1 + l‘^[{dx(t)Y + ^^ 0 ^] 


= 0 . 


(C.l) 


This equation admits an analytic solution similar to the one of the one dimensional 
case, namely 


y{x) 

y{x) 


isn {imy/Pci — lx + imy/Pci — Ic 
Im 

zsn (m\/Pci — lx + m^/Pci — 

Im 


2 


J^\ 

1+1-1 ) 



I z+i A 

I Pci-l j 


(C.2) 
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This soliton solution is similar to the solution of the one dimensional case discussed 
above. It is to see from the equations of motion that the map I —)■ —and —)■ — 

maps the one dimensional eqnation to the solitonic two dimensional one. 

We next consider the truly two dimensional case, thought of as an approximation 
for the fonr dimensional system of colliding shock waves in the limit of zero width 
for the shock wave, and in a limit of azimuthal symmetry in plane of the shock. It 
is convenient in two dimensions to use light-cone coordinates = t ± x, with 

d+= d^+= + d^) d-= d^-= ^{dt-d^). (C.3) 

In these light cone coordinates 

d^(j)d^(j) = (0)^ — = 4(9+0(9_0 d^d^cj) = 0 — 0" = 4(9+(9_0. (C.4) 

We now dehne the following coordinates 


s = — x'^ = x~^x 


<1 = 


X 

x^ 


For these coordinates we hnd that 
d^4>d^4> = +4 s{ds(j)f - y (<9q0)^ 
and also 


df,df^(j) = -F 4 [< 9 ^(s< 9 ^ 0 ) - ^[dg{qdg 


(C.5) 


, (C.6) 


( 2 ( 0 ^ 0 )(( 9 ,, 0 )(( 9 ^( 9 '' 0 ) =)d^(t)d^{du(t)d''(j)) = l&[s{ds(t)) +‘2s\ds(t)) dl(t)+ 
^((9g0)2((9,0) - 2g2(a,0)(0g0)(0g0,0) + 2^((9g0)3 + 2^((9g0)2(0j0)].(C.7) 

Snbstituting (C.6) and (C.7) into the eqnation of motion (2.2) we get that for the 
variables s and q the equation of motion takes the form 

2 

4[(9^(s(9^0) - ^(9q0 - ^<9^0] + m^(p = 

4/W0 


- 8 / 


1 -|- l‘^[m?(fP‘ — 4[s((950)2 — ^{dgCpy]] 

4 ls (( 9 , 0)3 + 2 s^{ds(j)fd^(j) + '^{dg(j)f{ds(j)) 


1 -I- /^[m^02 — 4[s{ds4>y ~ ^(<990)^]] 
-2q‘^{ds(j)){dg(j)){dgds(j)) + 2^{dg(j)f + 2^{dg(j)f{d^(j)) 
1 -I- /^[m^02 — 4[s{ds(t)Y ~ ^('9<j0)^]] 


For the special case of 0(s, q) = 0(s) the eqnation of motion reduces to 
4(9s(-S(9s0) -I- m^0 = 4l^s{ds(t))" 


,2. _ ,,4/0 02 + 2[(<9.0) + 2S020] 


(C.8) 


1 -|- /^[m202 — 4[s((9s0)2]] ’ 
which is the same as (2.6), so it reduces to the equation of motion of Heisenberg. 


(C.9) 
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